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We present recent results from the UKQCD collaboration's dynamical QCD simulations. This data has fixed 
lattice spacing but varying dynamical quark mass. We concentrate on searching for an unquenching signal in the 
mesonic mass spectrum where we do not find a significant effect at the quark masses considered. 
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1. Introduction 

Computing resources available for simulations 
of lattice QCD are now powerful enough to in- 
vestigate unquenching effects. UKQCD has em- 
barked on a programme of studying these effects 
in the light hadron spectrum, static quark poten- 
tial, glueball spectrum and topological sectors. 
The purpose of this paper is to investigate un- 
quenching effects in the first of these areas. 

It is well known that the lattice cut-off is a 
function of both the gauge coupling, and dy- 
namical quark mass, m**^" (see e.g. Q). For this 
reason, the philosophy we have chosen is to sim- 
ulate at points along the "matched" trajectory 
in the {f3,m^'^°-) plane i.e. defined hyfixed lattice 
spacing, a. This then disentangles lattice spacing 
artefacts with unquenching effects. Any varia- 
tion of a physical quantity along this trajectory 
can sensibly be attributed to unquenching effects 
rather than lattice systematics. This work has 
been published in full in 



2. Simulation Details 

A non-perturbatively improved clover action 
was used with lattice parameters displayed in Ta- 
ble |l| with a volume of IB'^ x 32. The last four rows 
contain the parameters for the matched ensem- 
bles, while the top simulation explores a lighter 
quark mass. Further details of the simulation ap- 
pear in |g. 



3. Mesonic Sector 

We begin outlining our results by indicating 

1 • 71 runitary / n ^unitary ^ 

our quark masses via Mpg /My where 
My (PS) is the pseudoscalar (vector) meson mass 
and the superscript "unitary" refers to the pa- 
rameters m^'^"- = m"°' in Table |l|. As can be 
seen our quark masses are only modestly light 
(c.f. M]^™*'^'-!//^;^"^*'^'-!/ ^ 0.18 in nature). 

The results for the hyperfine splitting are 
shown in Fig. 1 where the experimental points 
are plotted assuming ro = (0.49 ± 0)fm. Note 
that there is a tendency for the lattice data to 
flatten towards the experimental points as the dy- 
namical quark mass decreases within the matched 
ensembles. However the "unmatched" simula- 
tion {k^°^ — 0.1355) shows an increased negative 
slope, presumably due to finite-volume effects. 
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Figure 1. The hyperfine splitting. 
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Table 1 

Lattice parameters together with measurements of a and quark masses. 
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#conf. 


a-ro [Fm] 


aj [Fm] 


j^unitary j -^unitary 


5.20 


0.1355 


2.0171 


208 


0.0972(8)+^ 


o.iwtt 


0.578t;3 


5.20 


0.1350 


2.0171 


150 


0.1031(09)tf 


0.115t| 


o.70ot;^ 


5.26 


0.1345 


1.9497 


101 


0.1041(12)+;; 


0.118t^ 


0.7831^ 


5.29 


0.1340 


1.9192 


101 


0.1018(10)if 




0.835t^ 



5.93 1.82 623 0.1040(03)1^ 0.1186;;^ 1 



The J parameter, defined as, 
dMv 



J = AL 



(1) 



was calculated using 3 approaches: (i) "Pseudo- 
Quenched" i.e. where k^'^^ is held fixed and the 
derivative in eq.(|l|) is w.r.t. variations in k™'; 
(ii) "Unitary Trajectory" i.e. where k^"^^ = 
and the derivative in eq.(|l|) is w.r.t. variations 
in both K™' and k^°^ combined; and (iii) "Chiral 
Extrapolation of Pseudo-Quenched" i.e. taking 
the J values from Approach (i) and performing 
the chiral extrapolation m*"^" — > 0. The results of 
all three methods are shown in Fig. 2, together 
with the experimental point. There is evidence 
that the lattice value for J approaches the ex- 
perimental point as the sea quark mass decreases 
(see approaches 1 and 3). Note also that the 
"unmatched" simulation [k^'^^ = 0.1355) has a 
smaller J value. This fact is related to the com- 
ment in the previous paragraph. 

We use two methods to extract the lattice spac- 
ing. The first is via the Sommer scale rp, and the 
second is directly from the meson spectrum as 
outlined in [Q. The values obtained from both 
methods are displayed in Table |l|. Note that aj 
is « 10% larger than a^. This could be due to 
the experimental estimate of ro = 0.49 fm being 
« 10% too small. 

4. Unquenching Effects in Meson Spec- 
trum 

In order to investigate unquenching effects in 
the meson spectrum, we define the quantity 
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Figure 2. The J parameter from the three meth- 
ods described in the text. 



where is the scale determined from the phys- 
ical quantity Mi. When Si^j = then the lat- 
tice prediction of with scale taken from Mj 
agrees with experiment. Thus (S is a good pa- 
rameter to study unquenching effects. We expect 
that 5i^j[l3,'m'"^°-) — 0{a?) since we are using a 
non-perturbatively improved clover action. 

In Fig. 3 we plot 5i^j for the matched datasets 
where j = tq, and i is the scale determined from 
the string tension, ct, and the two mass pairs 
(p, 7r) and {K,K*) following The x-axis 

in this plot is (aM™*"''^)-2 ^ l/m^^^ so the 
quenched data point lies on the ?/— axis. 

We see disappointingly that, in the case of the 
scale determinations from the meson spectrum, 
there is no significant tendency of S towards zero 
as m^^'^ 0, i.e. the quantity 5 does not give us 
an indication of unquenching effects. However, in 
the case of cr, there is a statistically significant 
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Figure 3. Si as described in the text. 

variation of (5 — *■ as m'^'^'^ 0. This imphes 
that we do see unquenching effects in the static 
quark potential. 

5. Chiral Extrapolations 

In we used 3 approaches to perform the 
chiral extrapolations of hadron masses: "Pseudo- 
Quenched" ; "Unitary Trajectory" ; and a "Com- 
bined Chiral Fit" . In this publication we will refer 
only to the last approach. We take 



= Ao + AiMps{K 



Bo + BiMpsin 
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using the nomenclature M = aM and the first 
argument of M{k^°^; k™') refers to the sea quark 
and the second to the valence quark. 

The results of these extrapolations are shown 
in Table We stress that this functional form for 
the extrapolation is not motivated by theory, but 
is used as a numerical analysis technique in order 
to test for evidence of unquenching effects. As can 
be seen from Table ||, the parameters Ai and Bi 
are compatible with zero (to 2(t) and therefore we 
conclude that there is no evidence of unquenching 
effects. 



hadron 
Vector 
Nucleon 
Delta 



^0 



Ai 
-0.00411 
0.00611 
-0.0021? 



Bp 
0.6111 
1.231^ 
0.911^ 



Bi 
.0151? 
-.00111 
.0211 



Table 2 

Fit parameters from the Chiral Extrapolations 



6. Conclusions 

This paper attempts to uncover unquenching 
effects in the dynamical lattice QCD simulations 
at a fixed (matched) lattice spacing (and vol- 
ume) and various dynamical quark masses. This 
approach allows a more controlled study of un- 
quenching effects without the possible entangle- 
ment of lattice and unquenching systematics. 

However, we see no significant sign of un- 
quenching effects in the meson spectrum. This is 
presumably since our dynamical quarks are rel- 
atively massive, and so the meson spectrum is 
dominated by the static quark potential. This 
potential is, by definition, matched amongst our 
ensembles at the hadronic length scale rg, and 
so any variation of the meson spectrum within 
our matched ensemble must surely be a "higher" 
order unquenching effect which is beyond our 
present statistics. 

We have however shown that unquenching ef- 
fects exist in the static quark potential, and other 
work (using the same ensembles) has shown in- 
teresting unquenching effects in the glueball and 
topological sector . 
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